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ABSTRACT 

The extreme variability event in the TeV emission of Mrk 421, recently reported by 
the Whipple team, imposes the tightest limits on the typical size of the TeV emitting 
regions in Active Galactic Nuclei (AGN). We examine the consequences that this 
imposes on the bulk Lorentz factor of the emitting plasma and on the radiation fields 
present in the central region of this Active Nucleus. No strong evidence is found for 
extreme Lorentz factors. However, energetics arguments suggest that any accretion in 
Mrk 421 has to take place at small rates, compatible with an advection-dominated 
regime. 

Key words: galaxies: active - galaxies : jets - BL Lacertae objects: individual: Mrk 
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1 INTRODUCTION 

The intense 7-ray emission from radio-loud AGN (and in 
particular from blazars) argues strongly, independent of 
other evidence, in favour of the role of relativistic beaming 
in the physics of this class of source, as proposed almost 20 
years ago by Blandford & Rees (1978). In fact, the require- 
ment that the 7-ray source is optically thin to the process 
of electron-positron (e*) pair production, sets strong limits 
on the comoving radiation density at frequencies above the 
pair production threshold. 

The (potential) target photons can be modelled within 
the 7-ray emitting region, as well as at all scales along 
the line of sight, from the very central core of these Active 
Galaxies to the diffuse background field. 

Here we concentrate on the physical constraints that 
these observations impose within the inner ~ pc (for con- 
siderations relative to the background radiation field see e.g. 
De Jager, Stecker & Salamon 1994). In particular, we discuss 
the implications of the recent observation of the large and 
unprecedented fast variability event in the TeV emission of 
the BL Lac object Mrk 421 (z=0.031) by the Whipple Ob- 
servatory (Gaidos et al. 1996). On 1996 May 15 the fiux 
above 350 GeV varied by a factor of 20-25 in about 20- 
30 minutes, with a doubling timescale < 15 min. Note that 
the whole fiare has been observed. This variation sets tight 
upper limits on the size of the TeV emitting region. 



2 PAIR OPACITY CONSTRAINTS 

The optical depth to e* pair production by photon-photon 
interaction, t^--/, for 7-ray photons of energy x (in uifits of 
the electron rest mass energy) is proportional to the com- 
pactness of target photons, cc L/R, where 7? is the length 
of the path of the 7-ray in the target photon field. More 
precisely 



y(x) ~ r]{a)aTn{x,)x,R 



(1) 



where 2;, is the energy (in units of rricC?) of target photons 
with number density n(a;,), ctt is the Thomson scattering 
cross section and 77(1) ~ 0.12, where hereafter an energy 
spectral index a ~ 1 is assumed (Svensson 1987). 



2.1 Radiation field internal to the 7— ray emitting 
region 

Let us consider the limits deriving from any radiation field 
present in the same region where the TeV photons are pro- 
duced. 

An upper limit to the region size r^, implied by the ob- 
served 7-ray variability, is i? = r.y < ctvar5(l + z)~^ , where 
5 = (F — — 1 cos 6^)""^ is the relativistic Doppler factor 
of the emitting plasma moving with a Lorentz factor V at 
an angle 6 with respect to the line of sight and tvar is the 
observed variability timescale We consider the target pho- 



t Note that the size r-^ is estimated from the 7— ray variability 
and does not need to assume co-spatial emission with variable 
* E-mail: celottiSsissa. it flux at other frequencies. 
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tons at the threshold energy for pair production, where the 
cross section is maximum. Taking into account the Doppler 
transformations of the observed 7-rays, the corresponding 
target photon frequency is Vir ~ 1-2 x 10^* '^tcV with 
utcv ~ 1.2 X 10'^'' Hz (corresponding to an energy of 0.5 
TeV). The condition of transparency for the TeV 7-rays, 
r^-y < 1, then translates into a limit on the radiation field 
at i/ii. By assuming that all the observed monochromatic 
flux -FobaCi^ir) (erg s~^ cm~'^ Hz~^) is indeed produced co- 
spatially with the 7— rays, the constraint imposes a lower 
limit on the relativistic Doppler amplification (e.g. Dondi & 
Ghisellini 1995; see also Gaidos et al. 1996), namely 

5>10Fl'(^iOC3 (2) 

where the observed values for tvar = 1000 tvar,3 s and 
z^ir-Fobs(fir) — 4 X lO"'^^ crg s~^ cm~^ have been substi- 
tuted. These correspond to the quiescent state reported by 
Macomb et al. (1995) (at ~ 1.4 x lO" Hz). 

A Doppler factor 5 ~ F > 10 would therefore be re- 
quired to account for the observed TeV emission if the 
electron/positron population emitting at TeV energies also 
emits the observed IR fiux. Simultaneous monitoring would 
be required to establish that (see the results of the 1994 
and 1995 multifrequency campaigns, Macomb et al. 1995, 
Buckley et al. 1996). 



2.2 External radiation fields 

Analogously, one can consider constraints both on and from 
any (quasi-isotropic) radiation field external to the 7-ray 
emitting region, whose photons can efficiently interact with 
the TeV radiation. 

To estimate the photon-photon optical depth in this 
case it is necessary to estimate the compactness of this pu- 
tative field, over a scale size which corresponds to the dis- 
tance (from the central compact object) at which the high 
energy emission takes place. A typical perturbation/shock 
would naturally form, following a change e.g. in the bulk 
Lorentz factor F of the plasma in the jet, at a distance which 
depends on the intensity and time profile of the velocity 
change. For a variation AF ~ F occurring over the typical 
timescale associated with the dimension Vs of the compact 
object where the jet presumably forms (rg is the black hole 
Schwarzschild radius), a shock would develop at a distance 

R = Rj ~ rsF^(l + z)~^ Q Note that there is a minimum 
characteristic timescale for AGN variability, which is given 
by the light crossing time of the compact object ~ Vs/c, and 
if we assume the observed 7-ray variability does indeed cor- 
respond to this minimum timescale, than the 7-ray emitting 
region would be located at Rj ~ c tvarF^(l -I- z)~^. In other 
words, if the base of the jet varies on ~ rg/c, so would the 
observed 7-rays, though the scale they come from is ~ F^rg. 



t Here for simplicity it is assumed that the Lorentz factors of the 
bulk flow and of the front of any perturbation coincide, as we 
do not expect that the two values of F could differ significantly 
- for the purpose of this work - unless the shock front moves at 
extremely relativistic speeds relative to the mean fiow. 



2.2.1 Infrared field 

If one simply assumes that a diffuse IR radiation field is 
present over scale with a total fiux of intensity equal 
to the observed -Fobs(t'ir), then the transparency condition 
imposes again a lower limit on the Lorentz factor (which 
determines R-^ and therefore the compactness of target pho- 
tons). This limit corresponds to an extremely high 

F> 1.2 X 10^1^ (^.iOC^s (3) 

which in turn implies that the TeV emission is produced at 
distances of the order of i?^. ~ 4 x lO^^Fobs(fir) t~^r,3 cm- 
While there is no obvious reason why the plasma can- 
not propagate with this extreme F, the above limit depends 
on the hypothesis that a photon fiux comparable with the 
observed one pervades the region. One can then reverse the 
argument, assuming a more 'traditional' value of the Lorentz 
factor, say F ~ 10 Fi (e.g. Ghisellini et al. 1993), and de- 
ducing an upper limit on the diffuse IR field. This requires 
that the external field, Foxt(i-'ir), 

(1^) i^cxt(/^iO < 7 X 10-Vobs(i/iO (4) 

which significantly constraints the existence and properties 
of an IR emitter, e.g. in the form of a dusty torus, in the sur- 
rounding of the central object ^. In terms of luminosity, this 
corresponds to !-'ir-£'cxt(j^ir) ^ 10*" erg s~^, if the IR emitting 
region is located at distances > R-y ~ 3 x lO^^tvar.sFi cm. 

2.2.2 Starlight 

Within the inner ~ pc of the active galaxy, starlight con- 
stitutes a further radiation field in the surroundings of the 
compact object. While it is not possible to exactly quantify 
its intensity, an order of magnitude estimate, based on the 
centre of our galaxy, suggests a starlight radiation density 
corresponding to a luminosity Lstar ~ 10''° erg s~^ within 
Rstar ~ 1 pc, with a spectrum peaking, say, at fewxlO'* 
Hz. This radiation however does not contribute significantly 
to the TeV 7-ray opacity, giving rise to an optical depth of 

only Tjj^lf ~ 10" Latar,40-Rs^ar,pc- 

2.2.3 Line emission 

Variable emission in weak broad lines have been detected in 
the spectra of some BL Lac objects (Stickel, Fried & Kiihr 
1993; Vermeulen et al. 1995; Robinson & Corbett 1996). In 
particular, Morganti, Ulrich & Tadhunter (1992) report the 
detection of a broad emission feature in Mrk 421, with a 
luminosity (in Ha and Nil) of Lhc+ku ~ 1-5 x lO**" erg s^^. 
While during the observation reported by Morganti et al. a 
broad Ha line is not clearly resolved, the detection of this 
component in the spectra of other BL Lacs would support 
the hypothesis that in Mrk 421 the reprocessed emission 
also is produced by high velocity gas, supposedly in a broad 
line region. By adopting, in analogy with radio-loud quasars 
(and in fact quite arbitrarily for BL Lacs), a covering factor 

§ In a recent paper Protheroe & Biermann (1997) examine in 
details the effect of pair opacity above a typical AGN torus. 
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of ~ 10 per cent and a Thomson optical depth tt of a scat- 
tering diffuse medium of the same order, one can estimate 
the corresponding (isotropic) intensity of a ionizing radia^ 
tion field localized within a typical J?blr- The compactness 
of this putative broad line emitting region translates again 
into a 7-ray optical depth T^'"5y ~ 2Lion,4i-RBLj^ ig- Also this 
limit is too weak to set strong constraints either on the ion- 
izing radiation field or the bulk speed of the emitting TeV 
plasma. 

It is worth noticing that the estimate of the pliotoioniz- 
ing flux responsible for an observed broad Ha line, requires 
a luminosity of the order of 10''° erg s~^, not far fi-om the 
limit derived in eq. (4). 

2.2.4 Disc emission 

There is little evidence of the presence of thermal emission 
contributing to the observed radiation in BL Lac objects, 
probably the strongest indication being given by the narrow 
(and possibly broad) line emission. In particular, the limit 
inferred in eq. (4), under the assumption that the plasma 
bulk Lorentz factor is F < 10, sets constraints on the emis- 
sion from any material accreting, in the form of a disc, onto 
the central black hole. These can be translated into limits 
on the fundamental parameters regulating the structure and 
radiative properties of an accretion disc, namely the black 
hole mass M, accretion rate M and viscosity. 

In view of the numerical complication involved in the 
self-consistent computation of the emission at IR frequencies 
for the full range of the parameter space of M, M, viscosity 
(e.g. Szuszkiowicz, Malkan & Abramowicz 1996 and refer- 
ences therein), we adopt an extremely simplified approach. 
The emission from the disc can be roughly described (in 
vL^v)) by a power law (with the standard spectral index 
~ 1.3), peaking at the energy roughly corresponding to the 
maximum disc temperature, namely between the optical- 
UV and the soft-X band and extending to the low energies 
where the emission from the external radii dominates. 

We therefore consider the limits in mass and accretion 
rate of a disc not exceeding the luminosity of Mrk 421 in the 
optical-soft X-ray band, ly^vLuv, and extrapolate the omis- 
sion from the peak frequency to the IR band as a power law 
of index 1.3. This extrapolated luminosity is then compared 
with the limit given by eq. (4). 

However it is also necessary to estimate the location 
(within the limits of the optically thick Keplerian a disc 
assumptions, Shakura & Sunyaev 1973), of the IR emit- 
ting region, Rir- In fact, that allows us: a) to determine 
self-consistently the solution of the disk equations; b) to 
check that the IR emitting region is within the outer self- 
gravitating radius; c) most important, to estimate its loca- 
tion with respect to the 7-ray emitting region. The last point 
in fact determines the interaction angles (and therefore en- 
ergies) of the target photons for the TcV 7-rays. In particu- 
lar, for Rii < R^i, we consider that the absorbed IR photons 
have been mostly scattered and isotropized by an external 
medium with Thomson optical depth tt ~ O.Itt.-i. This in 
turn implies that the IR disc luminosity, consistent with the 
observed TeV emission, can be up to ~ 10 times the limit 
derived in eq. (4). 

Let us express the mass m in solar units and the 
accretion rate in Eddington units, rh = M/M^ (where 



Me = Lec^^ and Le is the Eddington luminosity). If we 
require that this bolometric luminosity does not exceed the 
observed emission VuvL^Uuv), this implies rh rm ^ 8 x 10~^. 
By extrapolating from the peak of the energy distribution to 
IR frequencies, the parameters of a standard thin Keplerian 
disc are then constrained to rh < 5 x 10~^Trjr^j . In this 
limit we have taken into account that for typical m > 10® 
and rh < 1, the disc radius which dominates the IR emission 
is smaller that R~f (where a Lorentz factor F = 10 has been 
adopted). For this range of m, m, and a viscosity parameter 
> 10~^, the disc at Rir is not self-gravitating. 

Therefore unless either Mrk 421 harbours a black hole of 
mass < 10^ Mq or the Lorentz factor of the 7-ray omitting 
plasma is P ^ 10 or the TeV 7-ray emitting region is at 
much larger distances, we are forced to conclude that the 
accretion rate in this object has rh < 10~^ — 10~^. 

On the other hand, observations on VLBI scales lead to 
estimates of the power emitted in the form of jet kinetic lu- 
minosity exceeding ~ lO""" erg s^^ (e.g. Celotti, Padovani & 
Ghisellini 1997), suggesting that indeed Mrk 421 is harbor- 
ing a much higher mass object. This could be reconciled with 
the low radiative (quasi-thermal) power if any accreting disc 
is radiatively inefficient. The deduced limits on r'n arc then 
consistent with the accretion occurring in the advection- 
dominated regime. 

2.2.5 Conclusions 

The constraints derived in this section imply that if a region 
comparable in size to the TeV production site R^ is per- 
vaded by a radiation field of intensity comparable with the 
observed one, then an extremely high bulk Lorentz factor, 
F > 10^, is required in order to overcome the limits from 
photon-photon opacity. 

However, one can easily envisage alternative possibili- 
ties which would not set such strong limits on the relativistic 
motion of the plasma. Namely: 

a) any quasi-isotropic IR field could be produced on scales 
much larger than 7?.y. F ~ 10 would require Rir ~ lO^ity 

(~ lO^^ F?o ivar.3 Cm). 

b) any isotropic (non-beamed) component is much weaker 
than the observed luminosity, and a typical F ~ 10 requires 
that -Fcxt(j^ir) ^ 10~*-Fobs(^'ir). This result implies that any 
disc in the nucleus of Mrk 421 accretes at a rate rh ^ 10~^ 
for black hole masses M > IO^Mq. 

Finally, radiation of stellar origin and a possible photoioniz- 
ing continuum do not significantly contribute to the opacity 
for TeV 7-rays. 



3 ON THE TEV EMISSION MECHANISM 

Let us now briefly consider the implications of the observed 
variability on the radiation mechanism producing the TeV 
emission. It is widely (but not universally) believed that 
the high energy 7-ray component in the spectra of blazars 
originates from inverse Compton scattering of relativistic 
electrons/positrons on soft photons (e.g. Sikora 1994 for a 
review). The origin of this photon field however is still a 
matter of debate, it could comprise both synchrotron pho- 
tons (SSC) and an isotropic photon field external to the 
synchrotron emitting region. 
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The condition that, in the comoving frame, seed pho- 
tons of frequency i/o are efficiently scattered to TeV ener- 
gies in the Thomson regime, imphes that r7 ^ 2 x 10^ and 
i/o 2 X 10^^ Hz, where 'yme(? is the energy of the scattering 
electrons/positrons in the jet frame. 

The energy density of photons at Vo, as seen by 
the emitting plasma, can be estimated as U'^-^i ~ 3 x 
lO'^Lext.^TT.-ir^^f'j^^ 3 erg cm~^ where a typical V ~ 
lOFi has been assumed and the dimension corresponding 
to the distance of the 7-ray emitting region. From the ob- 
served luminosity. Lobs ~ 2 x lO'*'' Lobs, 44 erg s^^, the syn- 
chrotron photon energy density is of the order of ?7gyn ~ 
0.6 Lobs,44 tvar,3 ^rg cm"^. Therefore f/s'yn/C^ext ~ 

20 which implies that the scattering of external pho- 
tons can dominate the SSC emission only for 5 > 20. 

Finally, rough equilibrium of the total luminosities 
emitted by Mrk 421 in the synchrotron and inverse Compton 
components of the spectrum (peaking in the UV-soft X-ray 
and 7-ray bands, respectively) implies a similarity of the 
magnetic and radiation energy densities, if the two radia- 
tive components are emitted in the same spatial region and 
most of the inverse Compton scattering is in the Thomson 
regime. This requires typical magnetic fields of the order of 
a few Gauss. 



4 SUMMARY AND DISCUSSION 

Consequences of the extraordinary variability event detected 
by Whipple for the physical properties of the emitting 
plasma and the radiation fields which surround the com- 
pact object and the powering of the central black hole in 
Mrk 421 have been examined. 

We find that high (> 10) or even extreme (> 10^) bulk 
Lorentz factors are required for the source to be transparent 
to TeV 7-rays, unless the radiation fields within and around 
the 7-ray emitting region are much lower than inferred from 
the observed flux. 

It should be recalled that BL Lac objects in general 
show little evidence of strong diffuse/isotropic radiation, and 
therefore the above considerations are not (yet) evidence for 
the most extreme values of the Lorentz factor. However we 
note that due to the weak dependence of the Lorentz factor 
on the compactness in target photons, we do not expect 
that such high F can be easily derived, at present, by 7-ray 
variability measures. 

Still there are no reasons a priori to exclude such a pos- 
sibility. Clearly, the beaming angle corresponding to F as 
high as thousands would create serious statistical inconsis- 
tencies if these results would be extended to all sources and 
the effective opening angle of the jet was of the same order. 
Indeed it is likely that the jet plasma flows within an angle 
^ 10~'^F^^. Furthermore, it is reasonable to imagine that 
jets reach such speeds only in the very inner core and that, in 
Mrk 421, we are observing the jet virtually aligned with its 
axis. This possibility would be also consistent with the mild 
apparent superluminal velocity /3appC ~ 3.8 c observed in 
the radio components of this source (Zhang & Baath 1990). 

One should note that even the weaker of the constraints 
derived on F is significantly higher than typical Lorentz fac- 
tors in BL Lac objects, as derived from limits on the SSC 
flux, VLBI radio measurement and statistical arguments 



(e.g. Ghisellini et al. 1993). Furthermore for these F the 
TeV emission is likely to be mainly due to SSC. 

Finally, we briefly consider the global energetics. We 
have shown that the assumption of a more 'traditional' value 
of the bulk Lorentz factor implies such stringent limits on 
the radiation field surrounding the nucleus that any accre- 
tion is likely to occur at rates m < 10~^. Interestingly this 
is further supported by independent recent findings. In fact, 
following the ideas proposed by Rees et al. (1982) and Fabian 
& Rees (1995), it has been suggested (Reynolds et al. 1996a) 
that M87, the nearby FR I galaxy, and by extension FR I 
radio galaxies in general, host an active nucleus of very low 
radiative efficiency as a consequence of the accretion be- 
ing advection-dominated. Extrapolating this to the plausi- 
ble beamed counterparts of FR I galaxies, namely BL Lac 
objects (e.g. Urry & Padovani 1995 for a review), one conse- 
quently predicts that these sources also are powered by an 
advection-dominated flow which is responsible for their low 
radiative (disk) efficiency. The jet power could then be ex- 
tracted, analogously to M87, either from the spin energy of 
the black hole (Blandford & Znajek 1977) or as Poynting- 
dominated outfiow along the system rotation axis of an 
accretion-powered object (Blandford & Payne 1982). This 
power could be of the order of Lom ~ 10'*^(a/mg)^_B|M9 , 
where B2 oc rn_[^ is the equipartition field (e.g. Narayan & 
Yi 1995), rh = 10~^m_3, ac is the specific angular momen- 
tum and rUg = GM/c? the gravitational radius. |^ 

Indeed we note that BL Lacs, and in particular Mrk 421, 
are hosted in elliptical galaxies, consistently with high values 
of the black hole mass and the suggestion by the above au- 
thors that these systems could be accreting hot interstellar 
medium quasi-spherically at the (Bondi) accretion rate. 

This hypothesis could also be consistent with the lack of 
much reprocessing gas in their surrounding as well as a low 
ionization fiux (but note earlier comments on broad lines. 
Sec. 2.2.3). If BL Lacs correspond to the final evolutionary 
stage of sources accreting through a radiatively efficient geo- 
metrically thin disk, this would be in (qualitative) agreement 
with the redshift distribution of BL Lacs, which is biased to 
low redshifts. 

The jet power deduced above could be compared with 
that estimated from the density of emitting electrons and the 
inferred Lorentz factors on VLBI scales, which imply a jet 
power of about three orders of magnitude higher, strongly 
suggesting that an high mass black hole lies in the very 
central nucleus of Mrk 421. The high value of the kinetic 
jet power can however also be reconciled with luminosities 
~ 10*^ erg s"'^ if jets of BL Lacs (and FR I) are mainly 
composed of an e^ plasma, as recently suggested on other 
grounds (Celotti et al. 1997, Bodo et al. in preparation) and 
in particular for M87 (Reynolds et al. 1996b). 

As it is often the case, some of the possibilities dis- 
cussed above could be settled when more and higher quality 
data become available. In particular we emphasize the im- 
portance of observations of variability at very high energies, 
which lead to much stronger physical implications when per- 
formed simultaneously with other energy bands. 

^ Note that the modeling of the high energy spectral distribution 
of BL Lacs favors a rather low magnetic field (in the emitting 
region) compared to high power sources. 
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